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Abstract 

We provide three Fortran programs which evaluate the QCD analytic (holomorphic) cou¬ 
plings for complex or real squared momenta Q^. These couplings are holomor¬ 

phic analogs of the powers a{Q'^y of the underlying perturbative QCD (pQCD) coupling 
a{Q‘^) = as{Q‘^)/7T, in three analytic QCD models (anQCD): Fractional Analytic Pertur¬ 
bation Theory (FAPT), Two-delta analytic QCD (2(5anQCD), and Massive Perturbation 
Theory (MPT). The index u can be noninteger. The provided programs do basically the 
same job as the Mathematica package anQCD.m published by us previously, Ref. [1], but 
are now written in Fortran. 
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Program Summary 

Program titles: kYkPleyit .for, A2danQCDext.for, AMPText.for 

The main programs (^AFAPText .for, A2danQCDext .for, AMPText.for^ and the 
tar-gzipped file containing all three programs (^AanQCDextFOR.tar.gz^, available 
from the web page: 
gcvetic.usm.cl 

Computer: Any work-station or PC where Fortran 95/2003/2008 (gfortran) is 
running 

Operating system in which the program has been tested: Operating system Linux 
(Ubuntu and Scientific Linux), Windows (in all cases using gfortran) 

No. of bytes in distributed programs: 

103 kB (AFAPText. f or), 107 kB (A2danQCDext. f or), 236 kB (AMPText.for); 99 
kB (AanQCDextFOR.tar.gz) 

Distribution format: tar.gz 

Keywords: Analytic (holomorphic) QCD coupling. Fractional Analytic Perturbation 
Theory, Two-delta analytic QCD model. Massive Perturbation Theory, Perturbative 
QCD, Renormalization group evolution. 

Nature of problem: Calculation of values of the running analytic couplings 
Aiy{Q‘^] Nf) for general complex squared momenta in three analytic 

QCD models, where Nf) is the analytic (holomorphic) analog of the power 

Nf)/ny. Here, AuiQ"^', Nf) is a holomorphic function in the complex 
plane, with the exception of the negative semiaxis reflecting the 

analiticity properties of the spacelike renormalization invariant quantities D{Q^) in 
QCD. In contrast, the perturbative QCD power Nf)/ny has singularities 

even outside the negative semiaxis (Landau ghosts). The three considered models 
are: Analytic Perturbation theory (APT); Two-delta analytic QCD (25anQCD); 
Massive Perturbation Theory (MPT). We refer to Ref. [1] for more details and 
literature. 

Solution method: The Fortran programs for FAPT and 25anQCD models contain 
routines and functions needed to perform two-dimensional numerical integrations 
involving the spectral function, in order to evaluate couplings. In MPT 

model, one-dimensional numerical integration involving AfiQ^) is sufficient to eval¬ 
uate any A^iQ"^) coupling. 

Restrictions: For unphysical choices of the input parameters the results are mean¬ 
ingless. When is close to the cut region of the couplings real negative), the 
calculations can take more time and can have less precision. 
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Running time: For evaluation of a set of about 10 related couplings, the times vary 
in the range t ~ 10^-10^ s. MPT requires less time, t ~ 1-10^ s. 
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1. General remarks 


The included Fortran programs evaluate the holomorphic couplings AuiQ"^) in the dis¬ 
persive approach to QCD, where we consider three analytic versions: Fractional Analytic 
Perturbation Theory (FAPT) [2] (and references therein), 2S analytic QCD (2(5anQCD) 
[3] (cf. also [1, 4]), and Massive Perturbation Theory (MPT) [5] (cf. also [6, 7]). These 
programs are constructed for complex squared momenta in a similar way as 

the previously published Mathematica package anQCD.m, Ref. [1] (also on the web page: 
gcvetic.usm.cl).^ For additional details on the three analytic QCD models with a view to 
implementing them numerically in Mathematica, we refer to Ref. [1]. 

A complication appears in Fortran, though, because in the dispersion integrals for 
FAPT and 25anQCD couplings the polylogarithm function \jiu{z) for complex appears 
and it has not been implemented in Fortran yet (in Mathematica it is implemented as 
PolyLog[—z/, ; 2 ]). In the attached Fortran programs this function is calculated as the 
following integral [14]: 


Li—n— 


d 


n+1 


di 


( 1 ) 

where v = n + 5. For better stability, the numerical integration of this integral over ^ is 
implemented in the complex plane along a ray in the hrst quadrant. This then results in 
two-dimensional integration for the evaluation of FAPT and 2(5anQCD couplings 
This integration is performed numerically with vegas routine [15]. In this approach, 
overflow and/or underflow problems appear at the edges of the unit square of integration, 
which are dealt with carefully in the programs. In the evaluation of couplings A^iQ"^) 
in MPT model, one-dimensional integrations are needed because the integrand involves 
Ai{Q‘^/0 no functions 

The Fortran programs are self-contained, i.e., no additional packages are needed. The 
needed explanations and instructions on the input parameters, compiling commands, and 
the output form, are all given at the beginning of each Fortran program. The calculations 
are more time-consuming in FAPT when the complex squared momenta are close to 
the cut negative semiaxis. In all models (FAPT, 2(5anQCD, MPT), the number of active 
quark flavors A^j is a hxed input integer. 

We wish to point out that this program, in Mathematica form [1], was already applied 
without problems to the evaluation of the nonsinglet structure function F 2 in deep inelastic 
scattering in FAPT model [16]. 


^ The analogous Fortran programs, but only for real values of the squared momenta , were men¬ 
tioned in Ref. [8] (and were made available on the web page: gcvetic.usm.cl). The Fortran program for 
the couplings T„(Q^) for integer index n and real in APT [9] and massive APT (MAPT) [10] was 
provided in Ref. [11], based on the corresponding program in Maple [12]. The Mathematica program for 
evaluation of the general power analogs A^iQ^) in FAPT was provided in Ref. [13]. 
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2. Practical aspects of the program 

2.1. Input parameters 

The programs can get compiled with the simple “gfortran” command. For example, 
the program A2danQCDext.for can get compiled by writing in the directory where the 
program is stored: “gfortran -o A2dan A2danQCDext.for” and then executed with the 
command “./A 2 dan”. Before compiling, however, the input parameters should be typed 
into the program, at (two) places which start with the string “INPUT”. In all three 
programs the following input parameters need to be specihed: 

1. Nf (“Nf”), the number of active quarks; 

2. A^in (“Nin”) and (5in (“delin”) indices, where u = A^m + ^in is the index of the coupling 
Au] where Ajn = 0 , 1 , 2 ,3,4; and 0 < (5in < 1 . 

3. \Q‘^\ (“AbsQ 2 ”, in GeV^) and 0 (“ArgQ 2 ”, in radians), where = |(5^|exp(i0) is 
the (complex) squared momentum (Q^ = —q^)- 

In addition, in the programs AFAPText.for and AMPText.for, the scale AjVj. (“gL2MS”, 
in GeV) of the underlying pQGD coupling needs to be specihed. 

Further, in AMPText.for, the MPT squared effective mass of the gluon mgj (“gM2”, in 
GeV^) needs to be specified. This mass is of the order of magnitude ~ 1 GeV^ [5^7], 
and is assumed in our program AMPText.for to be constant at all momenta (and thus 
independent of and of Nf). Instead of the input scale Ajvj. in MPT, which basically 

determines the strength of the MPT coupling = OisiQ"^ + ^^gi;MS)/ 7 r, one 

might prefer to use the value of 7 rMi(M|) = as{M'^ + mgpMS). This then determines 
the value of the scale IS.Nf= 5 . The values of other scales (for Nf = 4,3,6) can then 
be obtained by applying the (3-loop) quark threshold relations [17] applied within the 
(analytic) MPT model 


= Ml 


M2^+M3 

6 


36 " 


-I-M4 


iL 

216 


131,2 


4 + ^ (-6793 + 281(iV; - 1)) + F 3 


( 2 ) 


Here, A) = Nf - 1) and M„ = Nf), and 4 = Inf/i^^/m^] = lnK 

where fUg = ntgirng) is the MS mass of the corresponding quark entering at the threshold 
squared momentum (k ~ 1). We recall that in MPT (and FAPT) we use for 

the underlying pQGD coupling the 4-loop MS running coupling. In Table 1 we present 
the values of the scales corresponding to various values of mgj and 7rMi(M|) in MPT. 
A similar Table for the values of the scales A^r^ was given for FAPT in Ref. [1] (Table 1 
there). 

In A2danQGD.for program for 2(5anQGD model, the scheme parameter C 2 (= ^//do) 
was set equal to the central preferred value C 2 = —4.9 (cf. Table 2 of Ref. [1]). 
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Table 1: The scales Apfj., written in MeV, in MPT analytic model, for various values of and = 

as{M^ + mg[;MS). The threshold parameter k = 2 was taken; in parentheses, the results with the 
threshold parameter k = 1 are given. 


ml [GeV^j 

TrA^^'^'\Ml) 

^6 

^5 

A 4 

A 3 

0.5 

0.118 

88.3 (88.4) 

208.6 (208.6) 

291.5 (291.0) 

339.8 (338.4) 

1.0 

0.118 

88.3 (88.4) 

208.6 (208.6) 

291.8 (291.3) 

343.1 (342.1) 

1.5 

0.118 

88.3 (88.4) 

208.7 (208.7) 

292.1 (291.7) 

346.5 (345.4) 

0.5 

0.120 

99.8 (99.8) 

232.9 (232.9) 

321.8 (321.2) 

371.1 (369.5) 

1.0 

0.120 

99.8 (99.8) 

233.0 (233.0) 

322.1 (321.6) 

374.8 (373.6) 

1.5 

0.120 

99.8 (99.9) 

233.0 (233.0) 

322.5 (321.9) 

378.5 (377.2) 

0.5 

0.122 

112.3 (112.4) 

259.1 (259.1) 

354.0 (353.3) 

404.2 (402.1) 

1.0 

0.122 

112.3 (112.4) 

259.1 (259.1) 

354.4 (353.8) 

408.0 (406.7) 

1.5 

0.122 

112.3 (112.4) 

259.1 (259.1) 

354.8 (354.2) 

412.1 (410.6) 


2.2. Output 

After the execution of the program, the results of each program are written in the 
output hie AFAPText.dat (or: A2danQCDext.dat, AMPText.dat). If the input index is 
u = N + S (iV = 0,1,2, 3 or 4; and 0 < 5 < 1), the output will consist of the following 
couplings: 

0\.N+5{Q'^), ^N+l+siQ'^), ■■■, -^ 4 + 5 (Q^), (3a) 

AN+siQ^m^LO): AN+siQ^m^^^LO), ...^AN+siQ^m^m- (3b) 


Here we recall that Au(Q‘^) are the logarithmic derivatives analytically extended to non¬ 
integer index n t-A- u, where for integer n these derivatives are 


A^m 


(3^-\n-l)\ \d\nQ^) 


(n = l, 2 ,...) . 


( 4 ) 


The couplings AuiQ‘^) are analytic (holomorphic) in Q^, and perturbatively we have 
Au ~ {as{Q^)/7iy. Further, the couplings A^iQ^ are the analytic (holomorphic) analogs 
of the powers (as(Qy/Try, and they are linear combination of the couplings Au+xiQy 
{K = 0,1,...). For example, AN+siQ^ at N^+^LO precision {N is integer; 0 < 5 < 1) 
is the following linear combination: 


M 

AN+siQ"^) = AN+siQ"^) + km{N + 5)AN+m+&{Q‘^)- (5) 

m=l 


The expressions for the general coefficients km{N + 5) were obtained in Ref. [18]. These 
coefficients involve the Euler T function T(z/) = {d/dyvy) and its derivatives T*'™^(z/). 
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In Fortran, and are calculated by the routine [19] which is used in our 

Fortran programs. 
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